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Abstract 

Aims. We investigate the structure of the Galactic Hii region G353. 2+0.9, by analyzing (sub-)mm molecular-line and -continuum 
observations. This region is excited by the massive open cluster Pismis-24. We study the detailed morphology, distribution, and physical 
parameters (column and volume densities, masses, temperatures and opacities) of the molecular gas and dust. We are also interested in 
the variation in these parameters across the photon-dominated region. 

Methods. We observed various molecules and transitions to derive the physical properties of the molecular gas through line ratios, and 
both LTE and non-LTE analyses. The physical properties of the gas were derived with a Bayesian approach for the non-LTE analysis. 
Based on the continuum data at 870 ^im, we derived the column density of molecular hydrogen from the surface brightness and thus 
molecular abundances from the molecular column densities. We determined the mass of the dust from the integrated flux. We also 
carried out the simplest possible analysis to identify the clump candidates for gravitational instability, determining their virial parameter 
a. 

Results. The total mass of the gas in the region is ~ 2000 M©, while that of the dust is ~ 21 M©. The presence of a velocity gradient in 
the region, with clumps with redder Vlsr nearer Pis-24 suggests that the expansion of the ionized gas is pushing the molecular gas 
away from the observer. We unambiguously identify the ionization front in G353. 2+0.9, at the location of which we detect an increase 
in gas density and temperature. Its location and position angle is consistent with Pis-24 being the main ionization source. Almost no 
molecular gas is found south of the ionization front, at the location of the intense, elongated continuum and atomic-line emission, 
strengthening the hypothesis that Pis-24 is associated with G353.2+0.9. We find at least 14 clumps at difi'erent positions and LSR 
velocities, and we determine their physical conditions. The typical excitation temperatures are in the range of about 10-25 K, while 
H2 column densities are in the range ~ 10^^ - 10^^ cm~^. From the non-LTE analysis, we derive kinetic temperatures in the ranges 
1 1 - 45 K (CS) and 20 - 45 K (CN). The H2 number density is typically around ~ 10^ cm-^ from CS and few x 10^ cm-^ from CN, 
with maxima above 10^ cm~^. The abundances of the molecules observed are found to vary across the region, and appear to be higher 
in regions further away from the ionization front. 

Key words. ISM: Clouds - ISM: Individual: G353.2+0.9, NGC6357 - ISM: Molecules - ISM: Dust, extinction - ISM: photon-dominated 
region (PDR), Open clusters and associations: individual: Pismis-24. 



1. Introduction 

NGC 6357 is a complex of Hii regions and molecular clouds 
that form a very active star-forming region in the Sagittarius 
spiral arm. Optical, radio, and infrared (IR) images of NGC 
6357 confirm that it contains a number of distinct Hii regions in 
differ ent stages of evolution (e.g., Fellietal. 1990; Massi et al. 
|1997| ). Figure [T] shows a large cavity or a collection of smaller, 
connected cavities in the region, delineated by ionized gas. Weak 
and diffuse Ha emission permeates this feature. G353.2+0.9 is 
the bright emission region north of the cavity, seen in Fig.[T] Just 
55'' south of G353.2+0.9, lies the massive open cluster Pismis 
24 (hereafter Pis-24; Pi§mi§|1959 ). This clust er is thought to be 
the main source of io nization of G353.2+0.9 ( |Massi et al.|1997| 
[Bohigas et a l.'2004). It contains at least ~ 20 early-type (OB) 
stars, plus 24 0-type candidates ( Wang et al. 2007 ) and includes 
three stars that are amongst the brightest and bluest known in the 
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Galaxy, of spectral types 03 5 III(f*), 03.5 If*, and 04 III(f+) 
| Maiz Apellaniz et ar]|2007| ). [Massey et"ar| ( |2001| ) derived a 
distance of 2.56 kpc and an age of ~ 1 Myr for this cluster, and 
assuming that the molecular material is associated with Pis-24, 
we consider NGC6357 to be at the same distance. The large cavity 
is unlikely to have been formed by Pis-24, because of its clearly 
oflf-centre position. The morphology and the size of the cavity 
seem to suggest that it was shaped by the winds and/or supernova 
events of one or more clusters ( Wang et al. 20 07 ) situated in the 
proximity of the centres of the smaller bubble-like structures. 



Massi et al. ( 1997) performed a detailed study of the molecu- 
lar emission associated with two of the Hii regions in NGC6357, 
to wit: G353.2+0.6, and G353.2+0.9. The latter region is the 
younger one and exhibits signs of the presence of recently formed 
massive stars [e.g., ultra-compact Hii regions (UCHii), embedded 
sources with infrared (IR)-excess]. 

Our study is focused on G353.2+0.9.|Felli et slV (1990) ob- 
served it with the VLA at /I = 6 cm with an HPBW of T5 
(Fig. 7a of |Felli et al.|1990l ). As for the Ha emission, the high- 
resolution interferometric radio continuum observations reveal 
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Figure 1. The 8 jLtm emission image of NGC 6357, taken 
from the GLIMPSE survey (http : //www . astro . wise . edu/ 
^irtf/, Benjamin et al.|2QQ3t ). The cavity is clearly visible at 
the centre of the image. G353.2+0.9 is the bright region at its 
northern border. The location of Pis-24 is also shown in the figure. 
The coordinates are referred to the epoch J2000. 



a very complex structure of the ionized gas, with a well-defined 
sharp boundary running east-west (the "Bar", in Figure|2]). The 
emission is characterized by a strong intensity gradient to the 
so uth, while showing a more gentle de crease to the north (Fig. 8 
in lFelli et al. 1 19901 ). |Felli et aLl ( [T99Ql ) found three UCHn (A, B, 
andCinFig.[2|). 

The band image in Fig. [2] shows that in the central part of 
the nebula there is an elephant trunk-like region of obscuration 
(clearly visible also in HST images; flester & Desch|2Q05j , with 
an UCHii region and IR source at its apex. This source shows a 
near-IR excess and X-ray emission: it was identified from HST 
observations to be in the evaporating gaseous globule (EGG) 
evolutionary phase ( [Hester et al.|1996 ), making it the first X-ray 
emitting EGG. This embedded object was classified as having a 
spectral type B0-B2 ( Wang et al. 2007). The elephant trunk points 
toward Pis-24 and is thought to be formed by the radiation and 
stellar winds from the OB stars in this cluster. The IR emission is 
brightest along the sides of the trunk and on the south-western 
side of G353.2-F0.9, facing Pis-24. 

mapped G353.2-h0.9 in CO(l-O) and 



Massi et al. 



( .1997) 



^^CO(l-O). These data were complemented with observations 
of other molecules and transitions along strips in the north- south 
direction, to determine variations in physical parameters across 
the photon-dominated region (PDR). [Massi et al.| ( |1997| ) found 
that G353. 2-^0.9 is a face-on, blister-type Hii region, with most of 
the molecular material behind the Hn region and to the north of 
it. lFelH et al.| ( |199 0) suggested that G353.2-h0.9 is not associated 
with Pis-24, arguing that the southern sharp boundary (the "Bar") 
is produced by local ionization caused by embedded sources, that 
are also responsible for the radio emission of G353.2-h0.9. They 
also concluded that the nebula is ionization-bounded to the south, 
implying that there are considerable quantities of molecular gas 
in the region south of the ionization front. Molecular-line obser- 
vations do not however support this claim: Massi et al.| ([1997[) 
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Figure 2. i^s-band image of G353.2-h0.9 obtained with SofI 
(Massi et al, in prep.). The actual ionization front is indicated by 
"IF" (see text). The locations of the "Bar", Pis-24, and the elephant 
trunk are indicated. The red and green squares mark the location 
of some early- type stars in Pi s-24 and the three UCHn regions 
identified by ( Felli et al.|1990 ), respectively. The coordinates are 
referred to the epoch J2000. 



found very little molecular emission at the location of the "Bar". 



Bohigas et al. ( 2004 |) found that this elongated structure has 
to have a considerable extent along the line-of-sight (1 - 5 pc, 
Bohigas et al.|2004 l). While this dimension is comparable to the 
extent in R.A., it is much larger than the extent in DEC. This 
suggests that the "Bar" is a layer of ionized matter seen edge- 
on. It could be caused by the interaction of the photoionized 
photoevaporative fl ow with the free wind of the Pis-24 stars 
( Healy et"ar][2004 |. This implies that the molecular gas in the 
region has already been swept by the stellar winds, i.e. the region 
south of the "Bar" should be nearly devoid of molecular material. 

The pres ent study follows up on the work described in Massi[ 
et al. (1997]), which constitutes a first step in the study of the 
interface between the Hii region and the molecular cloud. Our aim 
is to clarify the morphology of the region, by observing optically 
thin molecular lines (e.g. C^^O), and to confirm the absence of 
molecular material south of the main ionization front, which is 
identified as IF in Fig. |2] (see Sect. |3. 11) . This would strengthen 
the hypothesis of an association between G353.2-h0.9 and Pis- 
24. The physical conditions of the gas are derived by means of 
a non-LTE analysis for those molecules with several observed 
transitions or with hyperfine structure, while for the remaining 
molecules we assumed LTE (local thermodynamic equilibrium). 
The observation of the continuum at 870 jam allows us to infer 
the dust column density and mass, and thus to determine those 
of the gas, by assuming a gas-to-dust ratio. With the H2 column 
densities determined in this way, we were able to calculate the 
abundance for the observed molecules. 



2. Observations and data reduction 

The molecular-line observations were carried out between Sept. 
1 and 9, 1999, with the 15-m Swedish-ESO Submillimeter 
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Telescope (SEST; program 63.1-0189). This radio telescope was 
operational on La Silla in the period 1987-2004. The telescope 
was equipped with SIS receivers and a high-resolution acousto- 
optical spectrometer with a total bandwidth of 86 MHz, and a 
frequency resolution of about 42 kHz. The spectrometer was split 
into two parts, and we observed with two receivers simultane- 
ously, one at a lower and the other at a higher frequency [e.g., 
C^^O(l-O) and C^^0(2-l)]. Detailed information about the molec- 
ular transitions observed are reported in Table [T] The columns 
indicate the molecules and transitions observed, their rest fre- 
quency, resolution in frequency and velocity, beam FWHM, main 
beam efficiency, spacing between points in the maps, and typical 
Tub rms noise per channel of the spectra, respectively. 

Most observations were made in frequency-switching mode, 
with a switch-interval in frequency sufficiently small to make the 
emission appear in both the signal and reference cycles, but large 
enough to avoid overlap between them. 

The SiO(5-4), CN(l-O), CN(2-1), and CH3CCH(6-5) lines 
were observed with position switching. The emission of the last 
two molecules exhibits hyperfine structure, making it necessary 
to observe in position- switch mode. The telescope pointing was 
checked every three hours on the nearby SiO maser source AH 
Sco, and was found to be accurate to within 5''. The same source 
was also used as the oflf-position. The calibration was made using 
the standard chopper- wheel method described in |Kutn er & Ulich| 
|l98l| ). Our molecular-line maps are centred ona = 17^24"^45^.6, 
S = -34° 1 V2(yn (J2000), coinciding with the "Bar". The angular 
extent of the observed region is about 5' x5' for CS(2-1) and (3-2), 
while it is ~ 3' X 3' for the other molecules and transitions, with 
spacing between the raster point listed in Table [T] (Col. 7). The 
line intensities in this paper are expressed in terms of the main 
beam temperature, defined as Tub = T^/t]mb- Data reduction and 
analysis for molecular- line data were performed with CLASS, 
part of the GILDAS (Grenoble Image and Line Data Analysis 
Softwar^ package. 

The morphology and the distribution of the molecular gas 
can be investigated in greater detail by decomposing the emis- 
sion profile into single Gaussian components at diff'erent Vlsr- 
Given the limited number of velocity components, we decided 
to decompose the emission profiles at every position by fitting 
diff'erent Gaussian curves, starting from optically thin transitions 
(e.g. C^^O). The characteristics of the components identified in 
this way were then used as a template for the decomposition of 
the emission profiles of the other molecules and transitions. 

CLASS off'ers the possibility to fit lines with hyperfine struc- 
ture, such as those of CN, by specifying the relative intensity of 
the hyperfine components in the case of optically thin emission 
and assuming that their ratios have their LTE values. CLASS uses 
the optical depth as a free parameter of the fit, and gives it as 
output of the procedure. 



We retrieved a map of G353.2-F0.9 at 870 ^tm (345 GHz), 
taken with APEX from the ATLASGAL survey ( Schuller et al. 
2009| ). The rms noise in the map is ~ 100 mJy/beam, determined 
in three regions free of emission around the Hii region. These data 
were analyzed with MOPSIC, the evolution of MOPSI (Map On- 
off' Pointing Sky dip Image), which was developed by R. Zylka 
(Obs. de Grenoble). 



Table 1. Molecular transitions observed. 
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AV 




Vmb 


Sp. 


rms 




(GHz) 


(kHz) 


(10~^km/s) 


(") 




(") 


(K) 


Ci^O(l-O) 


109.7822 


41.7 


11.4 


47 


0.70 


25 


0.11 


Ci^0(2-1) 


219.5603 


41.7 


5.70 


24 


0.50 


25 


0.20 


C^^ 8(2-1) 


96.4130 


41.7 


13.0 


54 


0.75 


25 


0.10 


CS(2-1) 


97.9810 


41.7 


12.8 


53 


0.75 


50 


0.10 


CS(3-2) 


146.9691 


41.7 


8.51 


35 


0.66 


50 


0.13 


CS(5-4) 


244.9356 


41.7 


5.11 


21 


0.50 


25 


0.25 


H2CO(2i,2-li.i) 


140.8395 


41.7 


8.88 


37 


0.66 


25 


0.30 


CN(l-O) 


113.4910 


41.7 


11.0 


46 


0.70 


25 


0.11 


CN(2-1) 


226.8748 


41.7 


5.51 


23 


0.50 


25 


0.23 


CH3CCH(6-5) 


102.5401 


41.7 


11.2 


51 


0.70 




0.17 


SiO(5-4) 


217.1049 


41.7 


5.76 


24 


0.50 


a) 


0.18 



http : //iram. fr/IRAMFR/GILDAS/ 



Note — Observed at selected positions only. 



3. Results and discussion 

3.1. Morphology 

Figure pi shows the maps of integrated line-emission J Tmb^^ 
towardsTj353.2+0.9. The molecular emission never extends sig- 
nificantly below = 0'\ confirming the lack of molecular 
material south of the "Bar" and ruling out the possibility that this 
feature is an ionization front proceeding southward. 

The molecular emission is concentrated between -10 and 
-hi km/s with, in some cases, strongly varying emission profiles 
between adjacent positions, as clearly visible in Fig. 4] In Fig. [5] it 
is possible to identify many different clumps. In higher- frequency 
transitions, such as CS(5-4) and CN(2-1), several clumps are 
resolved into two or more smaller clumps, or show an elongated 
appearance. 

We found that the emission can be separated into 14 clumps 
in six velocity ranges (see Figure[5|, i.e. -7.3 < Vlsr ^ -6.1 (A 
and P), -6.1 < Vlsr < -4.7 (B, C and L), -4.7 < Vlsr < -3.3 
(D and O), -3.3 < Vlsr < -1.9 (E and F), -1.9 < Vlsr < -0.4 
(G and H), and -0.4 < Vlsr < +0.8 (I,M an d N). The names of 
the clumps do not correspond to those used in |Massi et al. ( 1997| ), 
because of the higher spatial resolution in our present work and 
a different method of analysis (Gaussian decomposition versus 
integrated emission in velocity bins). 

Figure [5] shows that the various clumps have slightly different 
position for different molecules and transitions, which could be 
the result of different excitation conditions, optical depth, chemi- 
cal, or resolution effects. Furthermore, some clumps (such as P 
and O) are clearly visible only in high-density, high-frequency 
tracers. This is particularly clear in CS (5-4), which has the high- 
est resolution. The emission of these tracers shows the presence 
of multiple high-density, small cores within a larger clump. 

The clumps along the ionization front tend to have redder 
velocities than the others. This is especially true for low-density 
tracers (cf. Fig. |5]). Clumps N, H, and F, and D, B, and A ex- 
emplify this behaviour. This can be understood by taking into 
account the radiative and mechanical action of the stars of Pis-24: 
when neutral gas is exposed to the intense energetic radiation of 
an early-type star, it becomes rapidly ionized near the surface of 
the cloud. This gas is heated to T ~ 10000 K (a factor of ~ 100 
with respect to cold, neutral gas in the cloud), consequently caus- 
ing a comparable increase in pressure, thus leading to a rapid 
expansion. However, the expansion towards the neutral gas is 
stopped by the presence of the dense material of the cloud. In 
the opposite direction, the low-density ionized gas cannot halt 
the expansion of this overpressurized gas. The ionized material 



3 



A. Giannetti et al.: Molecular clouds under the influence of massive stars in the Galactic Hii region G353. 2+0.9 



CS(2-1) -10,+4 km/s 



CS(3-2) -10.+3 km/s 



CS(5-4) -10, km/s 




100 -100 -200 

(q) 

0^^^5(2-1) -10,+4 km/s 
1 1 I I I I I I I I I I I I I I I I I I 




1 



100 -100 -200 
(b) 

C^®0(1-0) -11, +3 km/s 
1 1 I I I I I I I I I I I I I I I I I I 




100 -100 -200 
(c) 

C^®0(2-1) -11. +4 km/s 
I I I I I I I I I I I I I I I I I I I 




100 -100 -200 
(d) 

CN(I-O) -12, + 10 km/s 



100 -100 -200 
(e) 

CN(2-1) -12, +4 km/s 
r 




100 -100 -200 
(f) 

H2C0(22,i-1ij) -9, +2. 5 km/s 
' I ' ■ ' ■ I ■ ■ ■ ■ I ■ ' ■ ■ I iU 15 




1 



10 



- 5 



- 



100 -100 -200 

(g) 



100 -100 -200 
(h) 



100 -100 -200 
(i) 



Figure 3. Maps of the integrated emission J Tmb^^ of the different molecular species and transitions, within the whole velocity 
range of emission. The first contour is the 3cr level. Molecule, transition, and integration limits are indicated above each map. The 
beam size is indicated by the filled circle. The observed positions are marked with a cross. The contours levels (in units of K km/s) 
are, respectively (lowest (step) highest): (a) 0.41 (1.0) 12.41, (b) 0.51 (1.0) 11.51, (c) 0.89 (2.0) 12.89, (d) 0.4 (0.4) 2.4, (e) 0.47 (0.7) 
7.47, (f) 0.81 (2.0) 16.81, (g) 0.98 (3.0) 27.98, (h) 1.25 (2.0) 19.25, (i) 0.92 (1.5) 15.92. Coordinates are offsets (arcsec) with respect 
to 17^24^45^.6 S = -?>A°\V2Q[!1 (J2000). 



moves predominantly away from the cloud, having an equal and 
opposite effect on the cloud. 

Therefore, the molecular emission indeed suggests that the 
ionized, over pressurized gas in G353.2+0.9, observed by Bohigas 
et al.| ( |20(j4| ), is expanding, thus pushing the molecular material 
away from the observer. 

3.2. Temperatures 

The first method used to derive the excitation temperatures uses 
the ratio of main beam temperatures of the (2-1) to (1-0) tran- 



sitions of C^^O. We first resampled the spectra to the same ve- 
locity resolution, and then convolved the C^^O(2-l) data with 
a Gaussian beam to match the spatial resolution of the (1-0) 
transition. 

We calculated the ratio Ri of rMB(2-l) to rMB(l-O) and de- 



termined Tex, assuming optically thin emission, from ( ,Levreault 
[T9881 ) 

where T^^^^i the excitation temperature of the C^^O(2-l) tran- 
sition. Assuming LTE, Tgx is the same for all levels. Since both 
transitions were observed simultaneously, the systematic calibra- 
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C^^0(2-1) -20, + 15 km/s; -0.2, + 9 K 
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Figure 4. TmbCC^^ 0(2-1)) spectral map. Each frame corresponds to a single pointing with an offset in a and S indicated in the figure. 
The X and y axes of each spectrum range from -20 to 15 km/s and from -0.2 to +9 K, respectively. 



tion uncertainties should not aff'ect the temperature ratios, which 
are instead influenced by uncertainties in the Gaussian decomposi- 
tion of the emission profile. Taking into account the uncertainties 
in the Gaussian fit performed with CLASS, we estimated that the 
uncertainties in Tub should not exceed 5-15%, where both lines 
are detected above 3cr. Assuming a typical uncertainty in Tub 
of 10%, the relative uncertainty in Tgx is ~ 14%. We found that 
Ri is always between ~ 2 and ~ 4, indicative of optically thin 
emission throughout the region. The Tq^(C^^O) (derived from 
Eq.[T]) is fairly uniform, typically between ~ 15 K and ~ 25 K 
for all clumps. The higher values are found along the ionization 
front, while for clump C, which is associated with the elephant 
trunk, we have Tex ~ 20 K. 

Another temper ature probe is methyl acetylene (CH3CCH) 
( Bergin et al.||1994| ). This molecule has a high critical density, 
and its emission comes from high-density regions. Our observa- 
tions are limited to four positions: (0'',50'0 (clump C), {15'\15'') 
(clump E), (-75'M25'0 (clump B), and (0'',-25'0. We detected 
CH3CCH only at (0'',50'0 and (75'', 75^0, while the detection 
at (-75'', 125'') is uncertain. Only at the first position were four 
components of the K-ladder visible, giving a reliable temperature 
estimate. The Boltzmann plot analysis (Fig. [6]) gave a Tex ~ 45 K 
for clump C, which is higher than that derived from C^^O. At the 
other positions, we detected just two components, leading to very 
uncertain temperature estimates (22 K, clump E; 45 K, clump B) 
that, however, are roughly consistent with those determined from 
CS (see Sect. |3.8| ). This higher temperature could be due to the 
presence of internal heating sources in clumps B and C. 

Although the excitation temperature is formally only a mea- 
sure of the relative population of the energy levels in a transition, 
and therefore diflfers for diflferent transitions, under the assump- 
tion of LTE, it provides a fair estimate of the kinetic temperature. 



3.3. Column densities 

Assuming LTE conditions, the total column density of the molec- 
ular gas can be obtained from observations of the J ^ J - 1 
transition by means of the expression ( |Zielinsky|1999 Kramer & 
|Winnewisser| 1 99l\ 



1 3h Z 

e ^^ex 

SttV J 



kTex 



[J(T,,)-J(TnG))rx 



(2) 



where J(T) = (/zyjj_i/^)(e^^J'J-i/^^ - 1)"\ t]c is the eflSciency with 
which the antenna couples to the source, jd is the dipole moment 
of the molecule, Z is the partition function, Tbg is the background 
temperature, and Tub is the main beam temperature. 

Equation ^ implies that there is a linear relation between 
the column density and the integrated line intensity at fixed Tex- 
This equation holds only in case of negligible optical depth (r). 
However, from the detection equation 



(3) 



one can derive a correction factor t/(1 - e"^), which makes it 
possible to calculate the column density while r < 2, with uncer- 
tainties less than 15% ( Kramer|1988) . In this exp ressi on, r is the 
optical depth at the line centre (see Sect. 3A_ and 3.8). 

In Table [2] we list the average value of the column density in- 
side the 3cr contour of integrated intensity, derived from diflferent 
molecules and transitions, for each distinct identified component. 
The derived excitation temperature is also listed, but we note that 
the column densities are computed assumi ng T^ ^ = 20 K for all 
molecules and a r derived according to Sect. 3A_ The assumption 
of a constant excitation temperature does not aflfect the results, 
because its range is small (15 - 25 K) and the variation in column 
density is only of the order of < 15 - 20%. 
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(b) CS(5-4) 

Figure 5. (a) Integrated emission J Tmb^v of C^^O(l-O) superimposed on the image. Each plot corresponds to a different Vlsr- 
The Vlsr of the clumps is indicated in the respective panel. The dashed (blue) lines indicate the observed area. The (red) contours 
show the integral under Gaussian components fitted to the line profiles. The beam of the transition is shown in the first panel. The 
first contour in each panel is the 3cr level. The integrated emission levels (in units of K km/s) are (lowest (step) highest): (-6.7 km/s) 
0.17 (0.3) 2.27; (-5.5 km/s, C) 0.19 (0.3) 1.99; (-5.5 km/s, B) 0.20 (0.5) 4.70; (-4.0 km/s) 0.18 (0.7) 6.48; (-2.5 km/s) 0.19 (0.3) 
2.89; (-1.0 km/s) 0.18 (0.2) 1.78; (+0.2 km/s) 0.15 (0.1) 0.65. (b) As (a), but for CS(5-4). The integrated emission levels (in units 
of K km/s) are (lowest (step) highest): (-6.7 km/s) 0.19 (0.3) 2.59; (-5.5 km/s, C) 0.20 (1.0) 10.20; (-5.5 km/s, B) 0.19 (0.5) 3.69; 
(-4.0 km/s) 0.19 (0.7) 4.39; (-2.5 km/s) 0.18 (1 .0)8.18; (-1 . km/ s) 0.17 (0.2) 1.37; (+0.2 km/s) 0.14 (0.1) 0.74. Note that names 
of the clumps do not correspond to those used in jMassi et al.| ( |1997) ). 
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Table 2. Vlsr, position, mean LTE column density, excitation temperature (C^^O), FWHM of the lines, and diameter of the clumps. 



Cl. (Vlsr) 
(km/s) 


Offset 


(l0l4xcm-2) 


(l0l4xcm-2) 


A^H2CO,21 
(l0l3xcm-2) 


(K) 


Line FWHM 

(km/s) 


FWHM size 
(pc) 


A (-6.8) 


/ yi c 1 or\\ 

(-45,180) 


9.71 


8.88 


1.09 


20 


1.3 


0.35 


B (-5.6) 


(-75,155) 


18.50 


18.36 


2.09 


20 


1.8 


0.55 


C (-5.5) 


(0,40) 


8.14 


10.10 


3.79 


21 


2.0 


0.35 


D (-4.4) 


(-70,120) 


23.60 


20.47 


2.35 


20 


1.8 


0.56 


E (-2.3) 


(80,60) 


12.00 


12.08 


6.21 


15 


1.8 


0.46 


F(-2.1) 


(-50,100) 


7.43 


5.54 




17 


1.2 


0.40 


G(-1.2) 


(50,75) 


10.50 


10.87 


3.32 


13 


1.3 


0.27 


H(-1.2) 


(-50,75) 


7.47 


6.34 


1.15 


15 


1.8 


0.31 


I (0.4) 


(0,95) 


2.81 


3.12 


1.14 


15 


1.4 


0.30 


L (-5.6) 


(-75,95) 


5.98 


3.07 




21 


1.2 


0.34 


M (-0.3) 


(-75,125) 


2.86 












N (-0.3) 


(-50,50) 


2.94 




1.11 








(-3.4) 


(50,70) 




8.07 


6.56 


17 


1.1 


0.37 


P (-5.9) 


(25,25)^1^ 






1.20 









Note — The column density is determined by averaging the emission inside the 3cr contour in J Tmb^i^ of each molecule in each clump, assuming Tex = 20 K. In the first column, the 
mean Vlsr of the clump is also indicated. The notation Najj means that the total column density of the molecule A is derived from the transition ij. The position is derived from 
CS(5-4). 
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Figure 6. Boltzmann plot obtained from the CH3CCH data at 
(0'',50'0. 



3.4. Opacities and visual extinctions 

We can derive the opacities of the transitions of C^^O, from the 
detection equation (jsj). The C^^O lines are usually optically thin 
(r ~ 10"^), although they reach r ~ 0.1 and r ~ 0.2, respectively 
for the transition (1-0) and (2-1) in clump E, in the elephant 
trunk and in some of the clumps aligned with the bright emission 
to the west of the trunk, i.e. the ionization front. Assuming the 
standard value of ~ 8 for the abundance ratio X(^^CO)/X(C^^O), 
one derives values of r for the ^^CO of the order of 0.8 and 1.6, 
respectively, for the transition (1-0) and (2-1) at the same posi- 
tions. This result co nfirms that the em ission of ^^CO is marginally 
thick, as argued by Massi et aL] ( |199 7 ). The higher optical depth 
indicates that the molecular material has accumulated in these 
regions. Furthermore, non-negligible opacities could explain the 
slightly lower Tex south of the clump C, with respect to those 
observed for the ionization front. 

An estimate of the visual extinction can be obtained fro m the 
H2 column densities from the expression (B ohlin et al . 1 1 97 8 ) 



5.34 X 10-^^Nh2 mag 



(4) 



Typical values are 5 - 10 mag and 15-20 mag for C^^O and 
H2CO, respectively. The maximum values of Ay that we found, 
estimated from H2CO, were those of clumps C and E (~ 50 mag). 



3.5. Dust 

From the 870 jam image, kindly provided by the ATLASGAL 
projecpl we derived dust masses and densities, following 



Deharveng et al. (2009) 



Mh = 



S870M.m^^ 



^870M.mB870M.m(7'd) 



(5) 



where Sgvo^m is the total flux density at 870 jam, D is the distance, 
B870M.m(^d) is the emission of a blackbody with temperature equal 
to at 870 |Ltm, and Ky = Ko(y/yo)^ is the dust opacity p er unit 
mass at the indicated frequency. Qssenkopf & Henning ( T994| ) 
recommend kq = 0.8 cm g"^ at 230.6 GHz, fror n which we 
derived A:870M.m ~ 1-8 cm^ g"\ assuming JS = 2 (Hildebrand 



fl983 ). We assumed that the dust temperature is 30 K throughout 
the region. However, the resulting masses and column densities 
do not depend strongly on temperature, for = 20 K and 
Td = 50 K the diff'erence is a factor of three. The dust mass 
resulting from the total 870 jam flux is ~ 21 M©. A:870M.m has an 
uncertainty of a factor of two, which implies an uncertainty in the 
dust mass of the same factor. Furthermore, the uncertainty due to 
the dust temperature is of the same order of magnitude. 

The surface brightness Fy at 87 |j,m also allows one to d erive 
the gas column density. Following [Deharveng et al.|(2009| ) 



rF870M.m 



2 . 3mH A:870 M-mB 870 M-m ( ^d )^Z?eam 



(6) 



F870M.m is the peak value of the surface brightness, Qbeam is the 
solid angle covered by the beam, mn is the mass of a hydrogen 
atom, and we assumed a gas-to-dust ratio y = 100. The contribu- 
tion of the free-free emission, extrapolated from the radio data 
of |Felli et al.| ( |199Q| ) assuming optically thin emission, is less 
than a few percent. The rms noise of the map (~ 100 mJy/beam) 
corresponds to a column density of 3.1 x 10^^ cm"^ for T = 20 K, 
1.9 X 10^1 cm-2 for T = 30 K, and 1.0 x 10^^ cm'^ for T = 50 K. 

We decomposed the emission into four bi-dimensional 
Gaussian components with MOPSIC, plus three regions of difl'use 
emission: Fig. [7] shows the condensations identified. 

The morphology of the emission at this wavelength is very 
similar to that of the integrated molecular emission and to that 



http : //www.mpifr-bonn.mpg .de/div/atlasgal/l 



|et aL|2009 



SchuUer 
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Table 3. Gas masses derived from dust emission. 



200 



100 



-100 



-200 



Figure 7. APEX image at 870 |Ltm. The scale is expressed in 
Jy/beam. The APEX beam is shown as a filled circle. The loca- 
tions of the main components identified are show^n in the map. 
Component 5 was not fitted with a Gaussian, but was considered 
as a region of diffuse emission. 



of the regions of obscuration visible at near-IR wavelengths. 
The maxima in the molecular-line emission do not always co- 
incide with those of the dust. In particular, clump E shows a 
large displacement. We suggest that this could be an effect of a 
non-negligible r. 

To associate the dust components with molecular cores, we 
superimposed the map of the 870 |Ltm emission on that of the 
molecular emission, as shown in Fig. [8] The masses of the com- 
ponents identified in the APEX image and their associated molec- 
ular clumps are listed in Table |3] The total mass indicated in this 
table is computed by integrating the total 870 jam flux. The H2 
column densities derived from the dust emission are on average 
in the range 3 - 7 x 10^^ cm"^, while the maximum values are 
cm (components 1, 2, and 4). 

The visual extinctions corresponding to the maximum column 
densities are typically 100 mag. The average visual extinction 
lies between 15 and 40 mag, depending on the component. The 
corresponding mean and maximum volume densities, calculated 
assuming spherical symmetry for the clumps, are -10"^ cm"^ and 



Component 


Mass 
(Mo) 


Associated molecular clumps 


1 


370 


CI 


2 


180 


C2 


3 


210 


H, F 


4 


340 


E 


5(1) 


360 


D, A,B 


total^2) 


2100 





Note — The gas-to-dust ratio was assumed to be 7 = 100. The total mass is computed by 
integrating the total 870 |am flux. Not a Gaussian fit, but a region that comprises all the 
emission at the corresponding location. computed integrating the whole 870// flux. 



Table 4. Single-clump and total gas masses. 



Cl. 


ci8o(i-o) 
(Mo) 


cl8o(2-l) 

(Mo) 


H2CO(2i2-ll 1) 
(Mo) 


A 


130 


160 


170 


B 


320 


310 


330 


C 


270 


270 


350 


D 


570 


400 


300 


E 


210 


240 


300 


F 


70 


50 


140 


G 


130 


90 


60 


H 


180 


170 


250 


I 


20 


30 


60 


L 


20 


40 




M 


< 10 






N 


20 




20 







50 


90 


P 






90 


Total 


1950 


1810 


2090 



Note — The mass is computed integrating the column density within the 3cr contour in 
/ Tmb^i^ for each clump. The transition from which we derived the mass is indicated in the 
top row. If the clump is divided into sub-clumps (indicated by numbers; cf. Fig.jSj the mass 
reported is the sum of all the sub-components. 



abundances vary in the region, in a similar way for the three 
molecules. The abundances are lower near the IF and the elephant 
trunk, while they appear to increase further away from IF, away 
from Pis-24. The C^^O abundance varies from ~ 9.0 x 10"^ in 
the region of the IF, to ~ 1.9 x 10"^ at the location of E, to 
~ 2.3 X 10"^ more to the north, at roughly the location of D. In 
the same three regions, the abundances of ^^CO and H2CO are 
~ 1.0 X 10-^ ~ 1.5 X 10-^ and ~ 2.2 x 10"^ and ~ 6.0 x 10"^^ 
~ 1.5 X 10"^, and ~ 1.1 x 10"^, respectively. We expect that 
this variation is at least partially due to an increase in Tk at the 
location of the ionization front. A different Tk can modify the 
abundance through the H2 column density, which was derived 
from the 870 |j,m emission. 



from H2CO, under the assumption of LTE. The low-density layers 
of the clumps are not visible in the 870 |^m map owing to its high 
rms noise. 

3.6. Molecular abundances 

In Fig.|9j we show the abundances of ] ^CO (not in the present 
dataset, taken from |Massi et al.||l997| , C^^O and H2CO as a 
function of position. To derive these maps, we smoothed the 
APEX image to match the resolution of the molecular transition 
considered and we divided pixel per pixel the molecular column 
density map by the H2 column density map, obtained from the 
APEX image, using Eq. |6] We set a threshold of 3cr in both 
molecular and H2 column densities. We can clearly see that the 



3. 7. LTE masses and volume densities 

The masses of the clumps are listed in Table |4] They were de- 
rived averaging the emission inside the 3cr contour to determine 
the average column density and then integrating over the beam- 
corrected area of the emission. We also took into account the 
optical depth r (determined from the detection equation) and a 
correction for helium, which contributes with a factor 1.36. The 
total mass of the complex is ~ 2000 M©, in agreement with 
that found from the 870 jam continuum emission. This mass can 
be compared t o the total mass of ionized hydrogen (290 M©, 
[Bohigas et al.|[2QQ 4) and to the total dust mass (21 M©, see 
Sect. 3.5|). The masses of single clumps range from a few tens to 
several hundreds of M©. 
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Figures. CS(5-4) emission (black and gray contours) superimposed on the 870 yum emission. The beam size of CS(5-4) is indicated 
by the filled circle, while the open circle indicates the APEX beam size. Each plot corresponds to a diff'erent Vlsr- The Vlsr of the 
clumps is indicated in the respective panel. The contours (in units of K km/s) for each clump are (lowest (step) highest): (-6.7 km/s) 
0.19 (0.3) 2.59; (-5.5 km/s, C) 0.20 (1.0) 10.20; (-5.5 km/s, B) 0.19 (0.5) 3.69; (-4.0 km/s) 0.19 (0.7) 4.39; (-2.5 km/s) 0.18 (1.0) 
8.18; (-1.0 km/s) 0.17 (0.2) 1.37; and (+0.2 km/s) 0.14 (0.1) 0.74. The lowest contour corresponds to the 3cr level in J TuBdv. 



The volume density n is determined assuming that the extent 
of the clump along the line-of- sight is equal to its size in the 
plane of the sky. The ratio of volume densities determined from 
low- and high-density tracers may be up to a factor of ~ 100, 
e.g., the mean volume densities derived from C^^O are typically 
~ 10^ cm-^ while those of H2CO are ^ 10^-10^ cm'^ The 
volume densities derived in this way are always of the same order 
of magnitude of the critical density of the molecular transition 
used to determine it. 

3.8. Non-LTE analysis 

The three transitions of CS, plus €^"^8(2-1), allowed a non-LTE 
analysis at ofi'sets where data points were available for all transi- 
tions. For this purpose, we used the statistical equilibrium, radia- 
tive transfer code RADE}^ for the approximation of a uniform 
sphere. The model requires as input the kinetic temperature, the 
H2 number density, the molecular column density, and the FWHM 
of the line. We assumed the Solar value of 22.5 for the ^^S/^^S 
isotopic ratio. The model returns the brightness temperature of 
the fines, the opacity, and the excitation temperature. The FWHM 
of the CS line was determined by averaging those of the various 
observed lines, at one position. The brightness temperature of the 
line was estimated using 

= ^;^^MB, (7) 



^ http://www. strw. le idenuniv.nl/$\sim$moldata/radex. | 
|html^|van der Tak et al.|2007l 



where §b is the beam size and §s is the source size. We assumed 
that all the transitions come from the same region, estimating 
the source size from the mean FWHM dimension of clumps in 
CS(5-4), that has the highest angular resolution. We varied the 
kinetic temperature between 10 K and 200 K, and the molecular 
hydrogen number density between 10^ cm"^ and 10^ cm"^, both 
in 50 equally spaced logarithmic steps. The column density was 
varied between 5.2 x 10^^ cm"^ and 5.5 x 10^^ cm"^ for CS and 
between 2.3 x 10^^ and 2.4 x 10^"^ cm'^ for C^'^S, in 44 equal 
logarithmic steps. 

To analyze the model results, we used a Bayesian approach. 
The Bayes theorem states that 



P(Tk, N, n\D, m) = -P(D\m, Tk, A^, n) P(Tj^, N, n), (8) 

where P(T}^,N,n\D,m) is the probability of the parameters 
T]^,N, and n, given the data and the model, called posterior, 
P(D\m, Tk, N, n) is the probability of the data given the model 
and its parameters, or likelihood, and P(T}<i, N, n) is the probabil- 
ity of the model parameters, which is called prior. The parameter 
^ is a normalization constant given by the sum of the individual 
probability of each model, in order to have the posterior normal- 
ized to 1. We used a constant prior for the model parameters, thus 
giving equal weights to every value of the model parameters. The 
probability of measuring a certain value for the intensity of a line 
is assumed to be represented by a Gaussian curve centred on the 
value obtained from RADEX for specific physical conditions of 
the gas, and with a a given by the uncertainty in the measured 
value. Therefore, we multiplied four Gaussian curves, one for 
each transition, and the probability density function (PDF) for 
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^^CO(I-O) 



C^^O(I-O) 






Figure 9. Maps of the abundance relative to H2 of ^^CO, C^^O, and H2CO, derived from i^CO(l-O) (top left), C^^O(l-O) (top 
right), C^^0(2-l) (bottom left), and H2CO(2i,2-li,i) (bottom right). The green contours shows the column density distribution of the 
molecules. 



Tc34s, computed with JAG;^ assuming that CS(2-1) is optically 
thick, to obtain P{D\m, T]^, N, n). The assumption of optical thick- 
ness for CS(2-1) is only used to derive the PDF for Tc34s, but is 
not used further in the analysis of the model results. Explicitly, 
the expression for P{D\m, Tk, N, n) is 



1 4 

^ ,--1 



(9) 



where the index / runs over the three CS lines and C^^S(2-1), // are 
the observed line intensities, yU/ are the modeled intensities, crij 
takes into account the rms of the spectrum and a 15% calibration 
uncertainty, and ^ is a normalization constant. 

To determine the physical conditions of the gas, i.e. the sin- 
gle parameters of the model, we had to integrate over the other 
parameters (marginalize) 



P(a\D, m) 



= J P(a,b,c\D, 



m) db dc. 



(10) 



From the PDF of the parameters, we derived the expectation 
values and the \cr range. The results for CS are summarized 
in Table [5] The columns shows the offset of the spectrum used, 
the clump name, the kinetic temperature and its 1 cr range, the 



http ://mcmc- jags . sourceforge .net/ 



number density of molecular hydrogen and its 1 cr range, the 
column density and its 1 cr range, the molecular abundance, the 
optical depth of CS(2-1), (3-2), (5-4) and C^^^{2-\\ and the 
excitation temperature of the same transitions, respectively. 

From these analyses, we found that the volume density of 
molecular hydrogen ranges between several x 10"^ and few x 
lO^cm"^, while Tk lies between ~ 11 K and ~ 45 K, for different 
clumps. The temperature derived from CH3CCH for clump C is 
higher than the kinetic temperature obtained from this analysis for 
(0'',50'0, while it is consistent with that at (0'',0'0, most probably 
because the CH3CCH beam takes in the emission from the region 
along the IF. However, CS(5-4) has its maxima at (0'',25'0 and 
at (25'',50'0 where we do not have data for CS(2-1) and (3-2), 
while at (0'',50'0 the emission of CS(5-4) is weak. 

There is a clear trend showing an increase in the densities 
towards the south, in the direction of Pis-24 and in the clumps 
aligned with the IF (Fig. |2]) west of the elephant trunk. This 
region shows typical densities of few x 10^ cm"^ and Tk ~ 
30 - 40 K. The temperature for the point at (0'',0'0 appears to 
be higher than the surrounding points, with Tk ~ 40 K, while 
(0'',50'0 has Tk ~ 18 K. The column density ranges between 
~ 2 X 10^^ cm"^ and ~ 2 x 10^^ cm"^, implying CS abundances 
of between 7.0 x 10"^^ and 4.0 x 10"^, derived from the ratio with 
the H2 column densities determined from the APEX data. The 
opacity of CS(2-1), CS(3-2) and CS(5-4) lies between 0.5 - 4, 
0.8 - 6, and 0.1 - 2.7, respectively. 



10 



A. Giannetti et al.: Molecular clouds under the influence of massive stars in the Galactic Hii region G353. 2+0.9 



The derived Tex-values are around 7-16 K for CS(2-1), 
6-11 K for CS(3-2), 5 - 9 K for CS(5-4) and 6 - 13 K for 
C^^S(2-1). The ratios of Tgx for different transitions are always 
within a factor of two of each other. 

A similar analysis was carried out for CN. In this case, 
P(D\m, Tk, N, n) was calculated by comparing total fluxes rather 
than line temperatures, due to hyperfine splitting, which is not 
included in the present model. Therefore we have 



P{D\m, Tk, N, n) 



1 2 



(Fi-Fm,i)2/(2£^2.)^Jj^-(Ttot,(l-0)-TIn,(l-0))/(2^T^) 



(11) 

where Fi is the measured flux, Fm,i is the output of the model, 
^tot,{i-^) is optical depth of the (1-0) transition as measured 
from the hyperfine satellite ratios, and Tmxi-O) is that predicted 
by the model. The uncertainty o-p^i takes into account the rms 
of the integral and sl 15% calibration uncertainty, and cTj is the 
uncertainty in Ztot, as given by CLASS. We used a Gaussian prior 
on Tk, centred on 35 K, with a (Xr^ = 30 K, given the results of 
the ana lysis carried out for CS and the CO results in |Massi et al. 
( 1997 ). The parameter (p is the normalization constant. Also in 
this case, the column and the number densities are constrained 
quite well, while the temperature is much more uncertain. 

The results obtained from RADEX for CN are summarized 
in Table [6] The columns shows the oflfset of the spectrum used, 
the clump name, the kinetic temperature and its 1 cr range, the 
number density of molecular hydrogen and its 1 cr range, the 
column density and its 1 cr range, the molecular abundance, the 
optical depth of CN(l-O) and (2-1), and the excitation temperature 
of the same transitions, respectively. The values listed in Table [6| 
for Tio are in very good agreement with those derived from the 
observations. 

Owing to the poor constraints on Tk, temperature maps are 
the most difficult to interpret. However, we obtain typical tem- 
peratures of between 25 K and 32 K for the whole region. The 
difference in found with CS at (0'^50'0 (18 K v^. 33 K) might 
be understo od if one considers that CN is a good tracer of PDRs 
( |Simon|1997 ), thus the emission may come from the more exter- 
nal material, directly heated and shocked by the interaction with 
the early-type stars of Pis-24. In all cases, these temperatures are 
usually within the Icr interval also found for E and at the edges 
of D and B, nearly facing Pis-24. 

In contrast, we derived a slightly lower temperature (~ 25 K) 
for C2 than for CI, which is nevertheless consistent with those 
cited above. 

The clumps in the region have a typical number density in 
the range ~ 1 - 6 x 10^ cm"^. The number density is highest 
along the IF (CI and C2), in D and in E, with values up to a 
few X 10^ cm"^. Clump E is particularly interesting, since it 
appears to have a resolved compressed layer facing Pis-24. The 
H2 number densities that we find are on average higher than 
those derived from CS, consistent with the idea that most of the 
emission comes from the high-density surface layers of the PDR. 

The column density of single clumps lies between ~ 5.0 x 
lO^^cm-^ and ~ 3.0 x 10i4^.j^-2 

After determining the column density for the various com- 
ponents, we were able to construct an integrated column density 
map, hence derive an abundance ma p similar to that of ^^CO, 

The emerging abundance 
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C^^O, and H2CO, as given in Fig. 

pattern closely resembles that of the other three molecules, indi- 
cating once more that the region of the elephant trunk and the 
brightest part of the IF is where the influence of Pis-24 is the 




Figure 10. Map of molecular abundance relative to H2, derived 
from CN data, through a non-LTE analysis. The green contours 
shows representative values of the column density, to clarify the 
distribution of the molecule. 

strongest. The abundance of CN is ~ 5.0 x 10"^ in the region 
of C, ~ 8.0 X 10"^ around D, and ~ 7.9 x 10"^ around E. We 
found that CN is enhanced around clump C and the IF, where 
the NIR emission is at its strongest and the ionizing radiation 
more intense: this is suggested by a less pronounced variation 
in abundance than that found from other molecules. For CN, the 
column density peak of the gas and that of the dust for clump E 
coincide, showing that its displacement for the other molecules 



can be caused by optical depth effects (cf. Sect. |3.6| or that C^^O 
and H2CO are frozen onto grains at the centre of the clump (given 
also the very low temperature derived from CS). However, a dis- 
placement is observed for the peak in the region of clumps A and 
B, again toward Pis-24. This and the temperatures comparable to 
those of CI might indicate that the IF/PDR extends here, even 
though the brightness at NIR wavelengths is much lower than 
nearer Pis-24. 

3.9. Virial masses 

The total mass of a spherical system in virial equilibrium is given 
by ( ,MacLaren et al.,1988j 

Mvir[Mo] - /:2^[pc]Az;[km/s]^ (12) 

where R is the radius and Av^ = 8 ln2cr^, assuming a Gaussian 
velocity profile and a density profile described by a power law 
p(r) oc R-'i with ^ < 3. The values of ^2 are 210, 190, and 126, 
respectively for ^ = 0, 1, 2 ( MacLaren et al.|1988) . The uncer- 
tainty caused by the unknown density profile is approximately a 
factor of two. This expression neglects the influence of magnetic 
fields, rotation, and internal energy sources, which are usually 
non-negligible in molecular clouds. 
To estimate R for Eq. 



12 



, we used the "eflfective radius", i.e. 
the radius of a circle with the same area as the clump above the 
FWHM level, corrected for the beam size. When the FWHM 
is smaller than the beam size, we assumed as an upper limit 
to the angular size, half of the actually observed FWHM. As a 
consequence of the large beam, some individual clumps may be 
blended and appear as one. For example, C is resolved into two 
different clumps in CS(5-4), and less clearly also in C^^O(2-l), 
while it appears to be unresolved in the other transitions. 

Figure [TT] shows the virial parameter a = Myir/MLTE as a 
function of Mlte, both determined from C^^0(2-l), where we as- 
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Table 5. Summary of RADEX results from CS data, for selected offsets. 



Offset CI. Tk Icr Icr A^cS [CS/H2] T21 T32 T54 T(^34s 21 ^ex ,21 T'ex ,32 T'ex ,54 ^exc345 21 

(") (K) (K) (lO^xcm-^) (lO^xcm-^) (lO^^xcm-^) (lO^^xcm-^) (lO"'^) ' (K) (K) (K) (K) ' 



(-100,200) 


A 


33 


10-64 


17 


4.8 - 36 


4.2 


3.9-4.9 


26 


0.7 


1.5 


0.5 


0.05 


14 


10 


7 


12 


(-50,150) 


B 


25 


11-35 


18 


8.3 - 30 


7.4 


6.1 -7.7 


5.2 


1.4 


2.5 


1.2 


0.09 


13 


10 


6 


11 


(0,0) 


C 


40 


10-72 


5.3 


1.9-8.3 


4.3 


3.1-6.1 


9.5 


1.2 


1.4 


0.1 


0.08 


8 


6 


6 


6 


(0,50) 


C 


18 


11-20 


25 


14-52 


21 


14-52 


2.7 


2.4 


4.1 


2.2 


0.17 


13 


11 


6 


10 


(-50,100) 


D 


26 


10-40 


22 


5.8-63 


16 


15 - 19 


3.4 


1.6 


3.1 


2.0 


0.12 


16 


12 


7 


13 


(100,50) 


E 


11 


10-15 


260 


130-480 


18 


15 - 19 


40 


3.8 


5.1 


2.7 


0.20 


11 


11 


9 


11 


(-100,100) 


F 


29 


10-40 


14 


5.8 - 25 


2.1 


1.6-2.5 


0.7 


0.5 


0.8 


0.1 


0.03 


11 


7 


6 


10 



Table 6. Summary of RADEX results from CN data, for selected offsets. 



Offset 


CI. 




lo- 


«H2 


\cr 


A^CN 


\cr 




T21 


,10 


T'ex ,21 


n 




(K) 


(K) 


(10^ X cm-3) 


(10^ xcm-3) 


(10^4 xcm-2) 


(10^4 xcm-2) 






(K) 


(K) 


(-50,150) 


A 


35 


21-57 


2.8 


1.2-7.6 


2.6 


2.2-3.4 


1.2 


4.0 


17 


10 


(-75,150) 


B 


25 


11-35 


2.3 


0.8-7.6 


0.7 


0.6 - 0.9 


0.9 


1.6 


9 


6 


(0,50) 


C 


33 


26-40 


18 


11-28 


2.3 


1.9-3.0 


0.2 


1.6 


67 


18 


(-50,100) 


D 


45 


30-68 


3.3 


1.7-6.3 


2.2 


1.6-3.0 


0.7 


3.2 


26 


11 


(75,50) 


E 


34 


28-42 


37 


23-69 


1.8 


1.4-2.2 


0.2 


1.1 


62 


25 


(50,75) 


G 


41 


28-58 


6.8 


3.0-16 


1.9 


1.4-2.5 


0.2 


2.2 


62 


14 


(-50,100) 


H 


33 


21-57 


2.5 


1.0-6.3 


1.7 


1.4-2.2 


1.2 


3.5 


15 


8 



a 




200 



Figure 11. Virial parameter a as a function of Mlte- ^vir and 
Mlte are both determined from C^^O. Mlte is calculated within 
the FWHM contour in J Tmb^v- The dashed line indicates a ~ 1, 
i.e. Mvir = Mlte. 



sumed that q = 2. The mass Mlte was derived within the FWHM 
contour in J Tmb^^- All clumps with masses above 50 M©, and 
also two with lower masses (M ~ 20 M© Gl and O) have a ^ I, 
thus indicating that these clumps might be gravitationally bound. 
We note, however, that this is an oversimplification of the problem 
of stability and must be taken with caution. 



3.10. Selective photodissociation 

G353. 2+0.9 is illuminated by several early-type stars. These stars 
deliver a huge quantity of energetic photons that dominate the 
chemistry and the heating of the gas in the PDR. The incident 
far-ultraviolet (FUV) flux is usually measured in terms of the 
ratio of the FUV flux in the region and the mean interstellar flux 
(the Habing Field, 1.6 x 10"^ erg cm"^ s"^). 

Considering the luminosities of the 03.5 If* [log(L/Lo) ~ 6], 
of the 03.5 III(f*) [log(L/L0) ~ 5.9], and of the 04 Ill(f-h) 
[logCL/Lo) ~ 5.8] stars given by |Weidner & Vink| (|2010), we 
can assume that these stars dominate the emission of energetic 
photons in the region. The fraction of luminosity emitted in the 
FUV band was estimated using a simple black-body law, between 
912 A and 2067 A. We obtained a total FUV luminosity of 



-FUV,tot 



and making use of 



= 1.2x10^ U 



1 



1.6 X 10-3 ergcm-2 s"! 47iD^' 



(13) 



(14) 



we found that Go ~ 5.6 x 10"^, for a representative projected 
distance of the stars from the elephant trunk of 0.66 pc. On the 
other hand, for the ionization front, we obtained Go ~ 2.0 x 10^, 
using a projected distance of 1.12 pc between the IF and the 
cluster. 

The high FUV flux may influence the ratio of C^^O and ^^CO, 
by means of selective photodissociation. This process tends to 
destroy the less abundant isotopologue, owing to the lower optical 
depth, thus less effective self- shielding capability. 

We obtained the ratio of the transitions of the two isotopo- 
logues (^^CO/C^^O) for both observed transitions, combining 
^^CO and C^^O values measured in positions less than half a 
beam apart. We corrected for the effects of opacity of the ^^CO, 
evaluated through the detection equation (given in Eq.|3]). The 
ratio of the Tub of the lines should therefore be representative of 
the relative abundance of the two molecules. The ^^CO/C^^O ra- 
tio is everywhere consistent with the standard value of the relative 
abundance for the two isotopologues (~ 8), with three exceptions. 
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The measured values of the ^^CO/C^^O ratio south of the ion- 
ization front, are a factor of two higher than the standard value. 
South of the elephant trunk, we did not detect C^^O. However, 
assuming emission from C^^O at rms intensity we obtain a lower 
limit to the ratio of 7 - 10. The model of |Visser QtH] ( [20091 ) 
shows that the relative increase in the ^^CO/C^^O ratio is about a 
factor 2-3, with the typical parameters of the region (Ay -5-10, 
Go ~ 10^ - 10^ and ~ 10"^ - 10^). The results reveal that selec- 
tive photodissociation does indeed occur south of the ionization 
front, and there is an indication that selective photodissociation 
is also taking place south of the elephant trunk. 

At the position of the elephant trunk, at the IF and for clump 
E the ratio is smaller than the standard value. In this case, the 
ratio may still be influenced by optical depth eflfects (in particular 
an underestimated opacity for ^^CO, which was determined from 
Eq.|3j assuming optically thin emission). 



3.11. Ionization front and geometry of the region 



Neithe r the present study nor that carried out by | Massi et"ar 
(1997]) found significant quantities of molecular material around 
the "Bar" (see Fig. [2] and [5]). This rules out the possibility that 
the "Bar" is an ionization front eroding a molecular cloud. This 
also implies that we cannot exclude an association of G353.2-h0.9 
with Pis-24 based on the position of this feature (as argued by 



Felliet al. 1990 prior to the availability of molecular line maps 



of the region). Energetic, spectral, and excitation analyses (e.g. 

1997; Bohig as et ari|2004 ) indicate that Pis-24 is 



Massi et al. 



indeed associated with G353. 2-^0.9 and that their proximity is not 
just a projection eflfect. The actual ionization front in G353.2-h0.9 
lies along the IR-bright feature labele d IF in Fig .[2| and is possibly 
associated with the UCHn region C ( |Felli et al.|f990) . 



We found a significant number of molecular cores along 
this bright ridge of emission and its continuation to the north- 
west, where the NIR brightness strongly diminishes. Furthermore, 
densities and temperatures are on average higher in this region. 
Our application of RADEX also revealed slight increases in Tk, 
volume density of H2, and opacity her e. There is intense radio 
continuum emission at 5 GHz (Felli et al.|199 0)) associated with 
this feature. This, together with the molecular gas distribution 
and physical conditions, confirm that IF is the main ionization 
front in G353.2-I-0.9. Here, the ionizing flux generated by the 
stars embedded in the region and by those of Pis-24, erodes the 
molecular cloud and pushes its material towards the north. When 
the resolution is high enough, one can see that the clumps near 
the ionization front, such as H, D, and C, are roughly parallel to 
it. High-density tracers are observed in rather small features at 
the edge of the ionization front, indicating the regions where the 
gas was compressed by the shock front. 

Molecular emission strongly decreases south of this ioniza- 
tion front and no massive clumps are observed here, although 
some emission i s visib le, especially in low-density tracers (see 
Fig. [5]). Bohigas et al. (2004 ) found that the region immediately 
surrounding the ionization front is characterized by a thin layer 
(10"^ pc) of very dense ionized material, where a photoevapora- 
tive flow is generated. 

We detect only very faint emission from CN(2-1) along the 
"Bar", which implies that there are very small clumps immersed 
in the PDR. 



3.12. Pismis-24 13(N36) 

This star is located in the northern part of G353. 2-^ 0.9 (cf. Fig. [2]) 
and is classified as a spectral type 06.5 V((f)) (Massey et al. 
[2001 ). Pismis-24 13 is worth noting because it seems to have 
produced its own Hn region in the molecular gas (see e.g.. Fig. 3 



m 



Hester & Desch|[2QQ5 ). This is confirmed by the radio con- 



tinuum and ion-line observations ( [Bohigas et al.|2004] ), which 



reveal free-free emission following very well the outer edge of 
the cavity and an increase in electron density in coincidence with 
this feature. This Hii region appears to be in the foreground with 
respect to the elephant trunk and ionization front. 

Dense star clusters often produce runaway OB st ars with high 
radial velocities through two diflferent mechanisms (Zinnec ker &| 
Yorke|2007| ): asymmetric supernova explosions and dynamical 
three-body en counters. The radia l velocities of these objects ex- 
ceed 40 km/s. [Gwamadze et al.1 ([2011 ) confirm that NGC 6357 
is rich in OB runaway stars ejected from the clusters within the 
cavity. Pismis-24 13 could be one of these OB runaways, ejected 
from Pis-24 in the direction of the molecular clouds, where its 
radiation and wind then created the observed Hii region. To con- 
firm this hypothesis we need a spectral measurement of the radial 
velocity of the star. 



4, Summary 

We have observed the Galactic Hn region G353. 2-^0.9 in several 
molecular lines, which has allowed us to distinguish at least 14 
clumps in its associated molecular cloud. We have determined 
temperatures, densities, and masses of each clump. We also identi- 
fied the location of the real ionization front in G353.2-h0.9. There 
is a tendency for the clumps near to the ionization front to have 
redder velocities than those further to the north. This is especially 
noticeable in low-density tracers (cf. Fig. [5]), and is caused by 
the expansion of the ionized, overpressurized gas pushing the 
molecular material away from the observer. 

Excitation temperatures derived from the ratios of the line 
temperatures of the C^^O(l-O) to the (2-1) lines indicate that Tex 
is in the range 15-25 K, with the higher values being found 
along the IF. The temperatures derived from CH3CCH, which is 
an eflfective tracer of kinetic temperature, are found to lie in the 
range of ~ 22 - 45 K. 

Assuming LTE and a constant Tex = 2 K, we derived the 
molecular column densities of ^^CO (from Massi et ar]|l997 ), 
C^^O, and H2CO, thus obtaining maps of molecular abundances, 
from their ratios with the H2 column density, which we derived 
from the APEX 870 jam image, assuming Tk = 30 K and a gas- 
to-dust ratio 7 = 100 (Fig.|9]). The maps show similar features, 
with a decrease in the molecular abundance in the region of the 
elephant trunk and the IF with respect to other region of intense 
molecular emission [~(75'',75'0, E; ~(-100'M25'0, D, B]. The 
molecular abundances derived in this way are uncertain by at 
least a factor of two, owing to the variations in Tk across the 
region. Nevertheless, the region of lower molecular abundance 
outlines the IF, clearly showing the area where the influence of 
early-type stars is the strongest. 

Column densities of molecular hydrogen derived from C^^O 
and H2CO, under the assumption of LTE and with the abundances 
calculated as described above, range between 10^^ - 10^^ cm"^. 
The visual extinctions are proportional to the H2 column density: 
typical values are in the range 5-30 mag depending on the 
transition used, while the maximum values, estimated from H2CO 
assuming LTE, are found in clumps C and E (~ 50 mag). The 
volume densities, determined assuming spherical geometry for 
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the clumps, lie between ~ 10^ cm" and ~ 10^ cm ^ . The volume 
density derived in this way is of the same order of magnitude as 
the critical density of the transition used to determine it. 

The total mass of gas in the region is ~ 2000 M©. Single 
clumps have masses in the range 10 - several x 10^ M©. The 
uncertainty in the mass for a given clump is dominated by the 
different physical conditions probed by the different transitions. 

A simple virial analysis shows that all the clumps with masses 
above 50 M©, in addition to two with lower masses (M ~ 20 M© 
Gl and O) have a ^ I, thus indicating that these clumps might 
be gravitationally bound. 

We performed a non-LTE analysis with RADEX, considering 
the four transitions of the two CS isotopologues in one case and 
the two transitions of CN in the other. 

To analyze the model results, we used a Bayesian approach, 
evaluating the likelihood P(D\m, Tk, N, n) based on Eqns. (|9| and 
( pTj ). We used a constant prior for all the parameters for CS and 
a Gaussian prior for Tk (yu = 35 K, cr = 30 K) for CN, according 
to our knowledge from CO, CH3CCH, and the results from CS, 
to reduce the degeneration on this parameter. 

For CS, we found Tk ~ 11-45 K, depending on the 
clump considered. The H2 number density typically ranges from 
several x 10^ cm"^ to few x lO^cm"^, but exceeds 10^ cm"^ for 
clump E, where Tk ~ 11 K. The CS column density lies be- 
tween ~ 2 X lO^^cm"^ and ~ 2 x lO^'^cm"^. Making use of this 
result, we determined the abundance of CS, taking the ratio of the 
molecular column density to that of H2 derived from the 870 jam 
emission. The abundances were found to lie between 7.0 x 10"^^ 
and 4.0 x 10"^ (Table|5]). These are lower limits because we do 
not know the column density of each velocity component, but just 
that of the strongest one. However, the results should not change 
by much, given that at a certain offset there is usually a single 
velocity component that dominates the emission. 

For CN, we were able to construct maps for Tk, N, and n. 
We found that Tk lies typically in the range of 25 - 32 K, with 
a maximum of ~ 45 K, but Tk is the most poorly constrained 
parameter, thus making the maps quite difficult to interpret. The 
H2 number densities that we found are on average higher than 
those suggested by CS, in the range of ~ 1 - 6 x lO^cm"^, which 
is consistent with the idea that most of the emission comes from 
high-density surface layers of the PDR ( |Simon|1997| ). The tem- 
perature Tk shows a similar behaviour: the temperatures found 
from CN are usually slightly higher than those derived from CS 
or C^^O. In the region of the elephant trunk and the IF, the vol- 
ume density is even higher, ~ few x lO^cm"^, which is similar 
to the values also reached for E and D. The column density lies 
between ~ 5.0 x lO^^cm"^ and ~ 3.0 x lO^'^cm"^. Having a map 
of CN column density for each clump, we summed them at every 
position, obtaining an integrated A^cn map, similar to those of 
C^^O, ^^CO, and H2CO. This allowed us to obtain an abundance 
map (Fig. [To]) that resembles those of the other molecules. The 
abundances found are in the range 5.0 - 7.9 X 10"^ (Table|6]). This 
map showed that CN is enhanced in the region of the trunk and 
along the IF, where the decrease in molecular abundance is less 
pronounced than for other molecules. 

We did not find significant quantities of molecular material 
in the region near the "Bar", which had been previously thought 
to be an ionization front. This is consistent with the idea that 
the "Bar" appears to be a layer of ionized matter seen edge- 
on (Bohigas et al.| |20Q4 ), the result of the free wind from the 
massive stars of Pis-24 interacting with the photoevaporative flow 
generated at the true ionization front (see Fig. [2]). The presence of 
a very faint feature in CN(2-1) along the "Bar" possibly suggests 



that the small quantity of molecular material in this region could 
be distributed in very small condensations inside the PDR. 

The high incident FUV flux strongly influences the shape and 
the properties of G353.2-h0.9. We investigated the presence of 
selective photodiss ociation of C^^O making use of ^^CO data of 



Massi et al. 



(1997). The ^^CO/C^^O ratio is nearly everywhere 



consistent with the standard value of the relative abundance of 
the two isotopologues. We found that south of the ionization front 
and the elephant trunk, C^^O is underabundant with respect to 
^^CO, being photodissociated by the strong, energetic radiation 
from Pis-24 stars. 

There seems to be a separate semispherical Hii region in the 
northern part of G353.2-h0.9, a ssociated with the s tar Pis-24 13 
(N36), which is an 06.5 V((f)) ( |Massey et al.|2001| . We propose 
that this star is a runaway O star, dissociating the molecular gas 
while making its way through it. Radial velocity measurements 
are needed to confirm this hypothesis. 
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